Impurity effects in unconventional density waves in the unitary limit 
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We investigate the effect of strong, nonmagnetic impurities on quasi-one-dimensional conventional 
and unconventional density waves (DW and UDW). The conventional case remains unaffected sim- 
ilarly to s-wave superconductors in the presence of weak, nonmagnetic impurities. The thermody- 
namic properties of UDW were found to be identical to those of a d-wave superconductor in the 
unitary limit. The real and imaginary part of the optical conductivity is determined for electric 
fields applied in the perpendicular directions. A new structure can be present corresponding to 
excitations from the bound state at the Fermi energy to the gap maximum in addition to the usual 
peak at 2A. In the dc limit, universal electric conductivity is found. 

PACS numbers: 75.30.Fv, 71.45.Lr, 78.30.Jw, 72.15.Eb 

I. INTRODUCTION 

The existence and behaviour of conventional (i.e. with constant gap) spin and charge density waves (SDW and 
CDW) is well documented^ . The thermodynamics of these systems was found to be very close to those of an s-wave 
BCS superconductor due to the similar, fully gapped density of states, but the transport properties are completely 
different. After the discovery of unconventional superconductors, the extension of the field of density waves (DW) 
into DW with wavevector dependent gap (termed unconventional) looks natural. In fact, after the earlier proposals 
in the context of the excitonic insulator^''^, this topic was rediscovered in the early 90's in various dimensions and 
systems^*^iSiL2i2ii£. Since then, the realization of unconventional or nodal density wavesii^i^ looks more and more 
possible: non-superconducting phase transitions without charge or spin ordering have been detected in a number of 
materials, and one of the possible explanations is provided by the unconventional density wave scenarioi^ii^*i^ii&. One 
of the main reasons of interest on UDW arises from high superconductors, where one of the competing models in 
the pseudogap phase is the d-density wave statepiii^. 

Recently we have studied the effect of impurities in the Born limit in unconventional density waves^^,. This treatment 
was justified from the fact, that this limit works very well for conventional density waves^*^, and the investigated 
physical quantities (for example the threshold electric field) showed convincing agreement with experimental data on 
a-(BEDT-TTF)2KHg(SCN)4^^i^^i^^. However, as is known from high Tc superconductorsSiS^, different impurities 
cause distinct effects on the same ground state: the Born and unitary scattering limit seems to describe Ni and 
Zn impurities, respectivel}i2&. From this, it looks natural to extend our earlier analysis on the thermodynamic and 
transport properties to the unitary limit. 

On the other hand, since conventional DW were mainly investigated in the Born scattering limit, it is instructive 
to study the effect of unitary scatterers on this state, partly to complete the picture, and partly due the interesting 
physics of this subject. 

In this paper we study impurity effects on quasi-one dimensional conventional and unconventional density waves 
at r = 0. The basic advantage of quasi-one dimensionality is that the nesting condition can be fulfilled at arbitrary 
fillings. First we examine the effect of resonant scatterers on conventional density waves. Interestingly, the density 
of states and the thermodynamics remains unchanged due to infinitely strong impurities, similarly to the effect of 
nonmagnetic impurities in s-wave superconductors in the Born limit^LSS. This surprising result follows from the fact 
that the nonmagnetic impurity enhances the rcnormalised order parameter A„ in the unitary limit like it docs in s- 
wave superconductor in the Born limit. As a result, a clean gap exists in the excitation spectrum for arbitrary impurity 
concentrations, and the low temperature physics is described by exponential functions with an activation energy. The 
unconventional situation gives more "conventional" results in the unitary limit. The thermodynamics looks very close 
to those of a d-wave superconductor in the unitary limit^"*'-^, and new localized states are visible around the Fermi 
energy. Similar phenomenon was observed in the density of states of isotropic p-wave superconductoriSSiffi, where a 
small island of states develops around the Fermi energy in the unitary limit. As a result of these new states at the 
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Fermi energy, depending on the direction of the appHed electric field and on the structure of the gap, new features 
are found for a; = A in the optical spectra along with the pair breaking peak at 2A, where A is the gap maximum. In 
general the gapless nature of optical excitations was detected experimentally in Q:-(BEDT-TTF)2KHg(SCN)4^, which 
coincide with our theoretical results, but for further conclusions more experiments are needed in the low temperature 
range. 



II. FORMALISM 



The single-particle electron thermal Green's function of DW is given hy^^^ 

G^^{k,iujn) = iujn - ^(k)/93 - /5icr3ReA(k) - p2CT3lmA(k), (1) 

where pi and a.i {i — 1, 2, 3) are the Pauli matrices acting on momentum and spin space, respectively, and for (U)CDW 
(T3 should be replaced by 1. A(k) = Ae*'^/(k), /(k) 1 in the conventional case and cos{bky) or sin(6fcj,) in the 
unconventional case. is the unrestricted phase (due to incommensurability) of the density wave, ^(k) is the kinetic 
energy spectrum, w„ is the fermionic Matsubara frequency. 

The Hamiltonian describing the interaction of the electrons with nonmagnetic impurities is given by 

^1 = ^ E e-"'^^*+(k + q)C/(R,)vI/.(k), (2) 



V, 



U{0) U{Q)e-'^^^ 



Rj is the position of the j-th impurity atom, Q is the nesting vector. 

The explicit wavevector dependence of the matrix elements^i*^ is neglected since no important changes are expected 
from it. Following the method of Ref. Il9l the self energy correction from impurities is given by 



ioJn) (4) 



where the the R index in Efi(ia-'„) means the position of an impurity over which the average should be taken, Ui is the 
impurity concentration. Here following the standard approach only noncrossing diagrams were taken into account. 
By fixing the ratio of C/(Q)/t/(0) and taking C/(0) to infinity, the self energy is given by 



T.{iujn) = -rii ( / ^-^G'(p,iw„) ) (5) 



We note here that in the special case of U{Q) — |[/(Q)|, the J7(R) matrix is singular, and the above calculations are 
not valid, but this condition corresponds to the fact that in real space the electron-impurity interaction is ultra short 
range, namely U{y) ~ <5(r), which is not the case in real systems. From Eq. jSJl, the self energy correction in the 
conventional case is obtained as 

SC^cO^-^^f (6) 



where g{Q) is the density of states per spin in the normal state at the Fermi energy, F = 2n,;/7rf;(0), u„ — tt>„/A„, a)„ 
and A„ are the renormalized frequency and gap: 

^« I 1- . ^ L (7) 



; + A2 

A = A„ I 1 - , ^ I . (8) 
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From this, the relation u„ = a;„/A„ = w„/A holds. In the unconventional case, the self energy correction is obtained 



as 



E(iw„)=r- ^— ^, (9) 



UnK -A 



where K{z) is the complete elliptic integral of the first kind. The gap remains unrenormalized due to the zero average 
of the gap function /(k) over the Fermi surface, m„ = a)„/A, and the Matsubara frequency is renormalized as 



a;„ = A«„-r- ^ ^. (10) 
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This is the same as in d-wave superconductors, the presence of backscattering (f7(Q)) drops out from the calculation 
and does not modify the result as in the Born limit. It is useful to introduce the quantity u„(a;„ = 0) = Co, which is 
determined from 



17Wtt) = 2a^^' (11) 



and will be used in further calculations. 



III. CONVENTIONAL DENSITY WAVE 

The density of states (DOS) is obtained as 

5(0) Vl - Vt^^ - A2 

where u — iun{iujn = oj + iS), 0(z) is the Heaviside function and the second equality follows from Eqs. lO and 
©. Hence the density of states remains unchanged in the presence of infinitely strong impurities, which is identical 
to the behaviour of s-wave superconductors in the presence of weak nonmagnetic impurities^LSS. As a result of the 
unchanged, gapped density of states, the thermodynamic properties such as the transition temperature or the specific 
heat remain the same as in the pure conventional density wave. This can be understood from the simple one impurity 
picture studied by Tiitto and Zawadowski in Rcf. Is^SSl The basic effect of impurities, the pinning comes from the 
interference between the Friedel oscillation and the density wave. For infinitely strong backscattering, however, the 
phase of outgoing electron is the opposite of the incoming one, hence the Friedel oscillation dies ou^. This simple 
picture has to be modified in the presence of the DW condensate, but the lack of Friedel oscillation still holds in the 
unitary limit, and no pinning is possible. 

In the transport properties, there are differences between the pure and impure systems. The optical conductivity 
for electric fields perpendicular to the chain direction still exhibits a clean gap for uj < 2A, but the divergent peak at 
2A turns into a sharp but finite cusp. 

IV. THERMODYNAMICS OF UNCONVENTIONAL DENSITY WAVES 

The density of states is obtained as 

^ ' = -Im-==K ^== = Im— -r-, (13) 



5(0) TT \VT~^J ^-A 

where u = iun{iujn = lo iS). It is identical to those of a d-wave superconductor in the presence of nonmagnetic 
impurities in the unitary limit^, and so does the thermodynamics as well, which can be borrowed from d-wave 
superconductors (Refs. I2^l25l and references therein). Consequently the change of the transition temperature is given 
by the Abrikosov-Gor'kov formula: 

-Infl^U^fi+pV^fiV (14) 



r/r, c^/A 



FIG. 1: In the left panel, A(0,r)/Aoo (dashed line), Tc/Tcq (solid line) and N{Q,V)/g{Q) (dashed-dotted line) are shown as a 
function of F/Fc. In the right panel, the density of states is shown for F/A = (dashed line), 0.01, 0.05, 0.1, 0.5 and 1 with 
increasing A''(0). The inset shows the localized state around the Fermi energy for F/A = (dashed line), 0.0001, 0.001, 0.005 
and 0.01 with increasing N{Q). 

where Tc and Tc^ are the transition temperature of the impure and clean system, respectively, p = T/2t:Tc, ^{z) is the 
digamma function. This formula holds also in the Born scattering limiliiS for both conventional and unconventional 
density waves as well as for unconventional superconductors in the presence of impurities considered either in Born 
or in resonant scattering limit^. The critical impurity scattering rate is obtained as 

27 4 ^ ' 

Using the parameters of Q;-(BEDT-TTF)2KHg(SCN)4, namely — lOK, vp — lO^m/s and lattice constant in the 
chain direction a ~ lO^^m, the critical concentration is estimated as rii ~ 0.001. In Fig. ^ we show the transition 
temperature, the residual density of states (i.e. N[Q)) and the zero temperature gap coefficient as a function of the 
scattering rate. The density of states exhibits localized state due to impurities around the Fermi energy superimposed 
on the usual gapless density of states of the pure system, which is manifested in the nonmonotonic nature of the DOS 
close to the Fermi energySS,, as it is shown in Fig. This state gives rise to a feature in the optical response as we will 
demonstrate it later. The identification of the localized or bound states is clearer for fully gapped systems: similar 
localized states were found in conventional density waves in the presence of one single, relatively strong impurity^^, 
and also the small island of states in unitary isotropic p-wave superconductors^Si^ around the Fermi energy signals 
the presence of new bound states. 

V. OPTICAL CONDUCTIVITY 

We calculate the optical conductivity for electric fields perpendicular to the conducting chain. In this case collective 
modes do not show up or can be neglected, depending on the explicit wavevector dependence of the gap^. Henceforth 
the optical response is calculated from the one bubble contribution, where self energy and vertex correction are taken 
into account in the noncrossing approximation. The real and imaginary part of the optical conductivity at T = is 
given bySSiS^: 

RcCTaa = — Re/(w), (16 

w Att 



oo 



ujlvt\aaa{io) = e^g{<S)vl^ [lvt\I{ijj) + 2 / lmF{u{x), u{x + uj))dx ] , (17) 
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where 



I{uj) — / {F{u{uj — x),u{~x)) — F{u{uj — x),u{—x)))dx 



(18) 



and Vx = vp, Vy = ^/2btb and Vz = V2ctc- In the following we discuss the different cases depending on the electric 
field orientation and on the gap. 
i. ^(k) = Acos{kyb), a — y: 



F{u,u') = — 

II' 



i E' { -uu' -l + ^\+K'{uu' 



2 



-Vr^(i?( 




(19) 



In the definition of the different F{ 
has to be used 



u,m') functions the argument of E and K is l/vT— while for E' and K' 

1/Vl ^ u'^ has to be used. In the present case, vertex corrections vanish similarly to the Born limit due to the 
mismatch of wavevector dependence of the velocity and the gap. In the real part a small peak develops close to w = 0, 
and moves to higher frequencies with increasing impurity concentration, but finally disappears as curves take more 
and more the form of a Lorentzian. Here the presence of bound states cannot be seen because the weight of scattering 
from the Fermi energy to the gap maximum is zero due to the zero velocity of quasiparticles at the latter point. In 
the imaginary part the cusp at a; = 2A smoothens as F increases, as it is seen in Fig. |5] The dc conductivity is 
calculated at T = as: 



dc,cos 2 



5(0)^;^ 



Att 



EJl + C^- 



7rF\ 
2A j 



(20) 



In the dc conductivities, the argument of E and K is 1/ y^l + C^. 





w/Aoo a;/Aoo 
FIG. 2: Real and imaginary part of the electric conductivity in the y direction for A(k) = Acos(6fcj;) 
of the reduced energy for different scattering amplitudes: F/A = (dashed line), 0.01, 0.02, 0.05, 
decreasing Recr(2A), lma{2A). 



are plotted as a function 
0.1, 0.2, 0.5 and 1 with 
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A(k) = Asin(fcj,&), a = y: 
1 



F{u,u') 



,'2 



I , 2 

-uu H \ 

2 I 3 3 



+ - 



T^:^/l-u'^^/l-u'^ 



1 



E'^Jl-u''^ - E^l-u^ 



-.K' - 



-.K 



2Auu'KK' {u + u') 



A2 



Tn 1 / \/l - 



2Au + m' \ w'ivT' 



(21) 



The real part of the conductivity exhibits a sharp peak at 2A and a small bump at A, indicating excitations from the 
localized state to the gap maximum for low concentrations. By increasing T, the former is suppressed and the latter 
becomes dominant. The imaginary part changes sign sharply at 2A, and a dip is present at A, as can be readily seen 
in Fig. O The dc conductivity is obtained as: 



vv 



''An^TTcl-AclE/r 



The conductivity is shown in Fig. |3| 




(22) 




0.5 1 1.5 2 2.5 3 3.5 



w/Aoo cj/Aoo 
FIG. 3: Real and imaginary part of the electric conductivity in the y direction for A(k) = Asin(bfcj;) are plotted as a function 
of the reduced energy for different scattering amplitudes: F/A = (dashed line), 0.01, 0.02, 0.05, 0.1, 0.2, 0.5 and 1 with 
decreasing Re(T(2A), increasing Im(j(A). 



iii. A(k) — Asm{kyb) or Acos(fcj,6), a — z: 



F{u,u) 



2(u'2 - u2) 



2^1 -u^E - 2 Vl - u''^E' + K' 



,u'(u — u') ^^u(u — u') 
' ] ' + K \ ' - 



(23) 



the vertex corrections vanish because the velocity depends on different perpendicular wavevector component (fc^) than 
the gap (ky). As F increases, the dominance of the A peak becomes more prominent than in the previous case in the 
real part of the conductivity. The imaginary part of the conductivity is zero for uj < 2A in the pure case, and exhibits 
a sharp peak at 2 A. The dc conductivity is obtained at T = as 



atl = 2e'g{0)v' 



E 



(24) 



w/Aoo w/Aoo 

FIG. 4: Real and imaginary part of the electric conductivity in the z direction are plotted as a function of the reduced energy 
for different scattering amplitudes: F/A = (dashed line), 0.01, 0.02, 0.05, 0.1, 0.2, 0.5 and 1 with decreasing Re(T(2A), 
increasing Imcr(A). 

The latter two cases seems to be consistent with experimental data on Q;-(BEDT-TTF)2KHg(SCN)4^ as far as the 

gapless nature of the optical response is considered, while the former with almost monotonically decreasing Recr(w) is 
different from the measured data. We refrain here from the evaluation of quasiparticle part of in chain conductivity, 
because the sliding collective mode associated with the phase of the condensate dominates this response^S,. We note 
however that the quasiparticle part of a-j^x(y>) is expected to behave very similarly to Uzz(y>)- 

The dc conductivities are shown in Fig. [S] at T = as a function of the impurity scattering parameter. In 
the perpendicular direction, the dc conductivities take the same value at the critical scattering parameter, namely 
e^(7(0)z;^ ^/Fc. Surprisingly, for small concentrations the dc response increases linearly with F, as opposed to the 
almost F independent behaviour in the Born limit This increasing behaviour is attributed to the fact that the 
creation of zero energy quasiparticlcs due to impurities is more efficient than the scattering of quasiparticles by 
impurities^. It is worth mentioning that in the dc conductivity the F ^ and a; ^ limit cannot be exchanged, as 
it is seen in Figs. |21 El El and |S1 this is why we obtain different dc conductivities in the pure case depending in the 
order of limits. However, we believe the right procedure is shown in Fig. El where the w ^ limit is taken first. The 
dc conductivity in all cases turns out to be universal^, since regardless to the scattering limit it takes the same value 
as F ^ 0, namely a^y^^"'"' = e2,g(0)w24/Aoo7r, a^l = g(}S)vl2 j I\qqt^ and cr^;;'"" = 0. The last equality holds since the 
electric current operator vanishes on the nodal points of the gap. 

VI. CONCLUSION 

We have studied the effect of nonmagnetic impurities in conventional and unconventional density waves in the 
unitary scattering limit in the standard noncrossing approximation. In the conventional case, no changes are found 
in the thermodynamics compared to the pure system, similarly to s-wave superconductors in the Born limit. In the 
presence of one single, infinitely strong impurity the Friedel oscillation disappear*^'^, since the phase of the incoming 
and outgoing electron is opposite. Consequently there is no interference between the density wave and the Friedel 
oscillation, and there is no pinning. In the presence of impurities with finite concentration, this simple picture seems 
to survive and the effect of impurities is canceled from the thermodynamics. 

As opposed to this, in the unconventional case, the thermodynamic properties are identical to those of a d-wave 
superconductor in the unitary limit. From the density of states, it is obvious that electrons are localized close to the 
Fermi energy, while at larger energies they remain almost unaffected by the presence of impurities (aside from the 
broadening of the lj = A peak). Also the change in the transition temperature is given by the Abrikosov-Gor'kov 
formula, which was also found to be valid in the Born limit Both the real and imaginary part of the optical 
conductivity seems to reflect the presence of localized states around the Fermi energy at certain gap structures by 
displaying a new bump at a; = A. This new feature seems to dominate over the pair breaking lo = 2 A peak as the 
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2.5 




r/r, 

FIG. 5: The dc conductivity is plotted at T = as a function of the reduced scattering rate for a cosinusoidal (sinusoidal) gap 
in the y direction: solid (dashed line) and in the z direction: dashed-dotted line. 



impurity concentration increases. We found universal electric conduction** in the dc limit. The comparison of the 
optical conductivity with experimental data seems to be difficult due to the lack of consistent investigations. This 
can be attributed to the fact, that the material which possesses most likely quasi-one dimensional UCDW ground 
state, the a-(BEDT-TTF)2KHg(SCN)4 salt enters into its new phase at lOK, and optical experiments below this 
temperature are very difficult. The only available data— reports some kind of pseudogap behaviour below Tc, which 
is compatible with our findings. Clearly, to make more decisive conclusions, further experiments are needed. 
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